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Quantum  coherence  of  a  narrow  band  particle  interacting 
with  phonons  and  static  disorder 

Qian  Niu  * 

Department  of  Physics 
University  of  California 
Santa  Barbara  CA  93106 


ABSTRACT:  This  work  deals  with  the  quantum  dynamics  of  a  narrow  band 
particle  interacting  with  phonons  cind  static  disorder.  The  present  theory  is  exact  in 
the  limit  of  small  band  width  compared  with  the  Debye  energy,  and  covers  various 
regimes  in  the  parameter  space  of  temperature  and  disorder  strength.  Therefore, 
the  theory  provides  a  unified  framework  for  studying  when  and  how  the  particle 
motion  changes  from  coherent  band  like  behavior  to  incoherent  hopping,  as  the 
temperature  and/or  disorder  strength  axe  increased.  The  theory  also  includes  the 
double  well  problem  as  a  special  case,  where  rather  complete  description  of  the 
particle  motion  is  obtained. 
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1.  Introduction 

In  this  paper,  we  study  the  quantum  dynamics  of  a  narrow  band  particle  mox'ing  in 
the  bulk  or  on  the  surface  of  a  crystalline  solid.  We  axe  interested  in  the  case  where  only 
the  lowest  level  of  each  potential  well  cam  be  occupied,  and  where  the  particle  moves  from 
one  well  to  the  next  through  quantum  tunneling.  We  aire  concerned  with  the  influence  of 
phonons  and  static  disorder  on  the  particle  motion.  In  particular,  we  want  to  know  when 
and  how  the  particle  changes  its  behavior  from  coherent  band  like  motion  to  incoherent 
hopping. 

This  is  an  old  problem  with  lasting  interest.  It  has  been  the  subject  of  small  po- 
laxon  theories  [1-7],  and  is  one  of  the  central  issues  on  quantum  diffusion  of  light  atoms  or 
muons[8-14].  In  the  general  context  of  quzintum  dynamics  with  dissipation,  the  present 
system  serves  as  a  fine  model  for  the  so  called  super-ohmic  coupling  category.  [15-21] 

The  main  objective  of  the  present  work  is  to  develop  a  general  theory  which  is  exact 
in  the  limit  of  small  band  width  compared  with  the  Debye  energy  of  the  phonons.  The 
theory  will  be  valid  for  essentially  the  whole  parameter  space  spanned  by  the  temperature 
and  disorder  strength.  It  therefore  not  only  unifies  earlier  theories  for  the  various  regimes 
of  the  parameter  space,  but  edso  make  it  possible  to  investigate  the  intermediate  regions 
where  transition  from  coherent  to  incoherent  motion  of  the  particle  takes  place. 

The  physics  of  transport  phenomena  usually  has  three  levels  of  description:  micro¬ 
scopic,  macroscopic  and  intermediate.  On  the  macroscopic  level,  one  observes  a  hydrody¬ 
namic  law  such  as  a  diffusion  equation  characterized  by  a  diffusion  constant  depending  on 
temperature,  isotope  mass,  etc.  On  the  microscopic  level,  one  tries  to  single  out  the  part 
of  the  microscopic  system  which  is  eventually  responsible  to  the  macroscopic  phenomenon 
and  which  is  describable  by  the  Schroedinger  equation.  These  two  worlds  are  connected  „ 
by  an  intermediate  level,  where  one  describes  the  system  in  terms  of  Markovian  equations 
involving  a  few  degrees  of  freedom.  Boltzmann  equations  and  master  equations  are  such 
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examples.  It  is  only  on  this  level  of  description  that  the  question  of  quantum  coherence 
is  releveint.  On  the  microscopic  level  everything  is  coherent.  On  the  macroscopic  level  ev¬ 
erything  is  incoherent  (excluding,  of  course,  particular  systems  such  cis  a  superfiuid).  We 
emphasize  that  quantum  coherence  is  not  equivalent  to  the  Markovian  nature  in  the  trans¬ 
port  equations.  It  has  something  to  do  with  the  quantum  uncertainty  between  position 
and  momentum.  In  this  paper,  the  pzirticle  motion  is  termed  'incoherent'  if  it  is  Maxkoviem 
in  position  basis.  The  usual  Boltzmann  equation  is  Mcirkoviaui  only  in  momentum  basis, 
and  therefore  describes  a  coherent  motion.  Further  refinement  of  this  notion  will  be  done 
in  the  text. 

The  derivation  of  transport  equations  usually  starts  from  the  microscopic  level,  follow¬ 
ing  totally  different  routes  for  cases  with  qualitatively  different  behaviors.  Unfortunately, 
such  an  approach  can  draw,  at  best,  scattered  pieces  of  a  whole  picture,  and  is  unable 
to  answer  questions  like  how  a  transition  from  coherence  to  incoherence  takes  place.  For 
example,  the  conventional  small  polaron  theory  follows  such  a  practice.  For  the  low  tem¬ 
perature  side,  one  starts  with  the  renormalized  Bloch  states,  and  uses  perturbation  theory 
to  determine  their  life  times.  When  the  inverse  life  times  are  smaller  than  the  renormalized 
beind  width,  the  particle  motion  is  said  to  be  coherent.  When  the  inverse  life  times  are  no 
longer  small  for  high  enough  temperatures,  it  is  said  that  the  original  picture  of  coherent 
motion  breziks  down.  The  theory  for  the  high  temperature  side  then  takes  a  totally  differ¬ 
ent  route:  one  starts  with  the  unperturbed  site  diagonal  states,  a'nd  then  uses  perturbation 
theory  to  calculate  the  rate  of  transitions. 

The  path  integral  approach  of  Leggett  et  al[19]  was  designed  to  give  a  consistent 
picture  of  the  particle  motion  in  various  regimes.  Their  non-interacting  blip  approximation 
works  remarkably  well  for  the  subohmic  and  ohmic  cases,  and  has  correctly  predicted  a 
biais  induced  transition  from  coherent  to  incoherent  tunneling  for  the  superohmic  case. 
The  approximation  fails,  however,  for  the  intermediate  region.  It  is  valid  in  both  the  zero 
bias  limit  and  the  large  bias  limit,  but  for  totally  different  reasons.  Recently,  Gorlich  et 
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al[20]  made  an  improvement  of  the  approximation,  giving  correct  result  also  for  moderate 
bias.  However,  it  is  still  not  clear  how  their  results  may  be  generalized  to  high  temperature 
cases  and  multi-site  problems. 

In  this  paper,  we  take  a  different  approach.  We  start  with  a  path  intcgred  represen¬ 
tation  of  the  transition  probability  for  the  particle  to  move  from  one  place  to  another  in 
time  t.  This  is  a  path  integral  defined  in  the  phase  space  of  the  particle,  i.e.,  the  space 
spanned  by  the  two  position  variables  appearing  in  the  particle  density  matrix[22].  The 
phase  space  path  integral  defines  a  phase  space  propagator  (or  Green  function).  A  Dyson 
equation  is  then  obtained  by  regrouping  the  terms  of  a  diagramatic  expansion  of  the  prop¬ 
agator  in  some  small  parameter  to  be  specified  below.  The  first  order  approximation  of  the 
self-energy  operator  then  reduces  the  original  many  body  problem  involving  the  phonons, 
essentially  to  a  single  particle  eigenvalue  problem  in  the  phase  space  of  the  particle[23-25]. 
This  provides  a  very  convenient  basis  for  further  explorations  of  various  regimes  and  their 
inter-connections. 

For  a  heavy  particle  such  as  the  hydrogen  atom,  the  tunneling  energy  splitting  is 
extremely  small  compared  with  the  Debye  energy  of  the  lattice  vibrations[12.13].  The 
ratio  of  the  two  energy  scales  provides  us  with  a  small  parameter.  Our  theory  is  based  on 
a  systematic  expansion  in  this  smeill  parameter.  Truncation  after  the  first  nontrivial  order 
should  give  a  very  good  approximation.  This  is  somewhat  in  line  with  the  philosophy  of 
the  small  polaron  theory  of  HoIstein[l]  which  also  treats  the  tunneling  energy  splitting  as  a 
small  parameter,  but  we  do  not  make  a  ■priori  assumptions  on  the  character  of  the  particle 
motion  (band  motion  or  hopping).  The  noninteracting-blip  approximation  of  Leggett  et 
al[19]  was  also  designed  to  take  advantage  of  this  smedl  parameter,  but  our  theory  is  more 
systematic. 

Equipped  with  our  general  theory,  we  have  studied  the  particle  motion  in  various 
regimes;  low  temperature  versus  high  temperature,  weak  static  disorder  versus  strong 
disorder,  and  two  wells  versus  an  array  of  wells.  We  have  rederived  the  transport  equations 
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in  various  extremes  of  the  parameter  space  and  revealed  once  again  the  fact;  the  particle 
motion  is  coherent  at  temperatures  low  compared  with  the  Debye  temperature  and  with  site 
energy  disorder  weaic  compared  with  the  renormalized  tunneling  energy  splitting;  whereas 
the  motion  is  incoherent  at  either  high  temperature  or  with  strong  disorder.  The  results 
on  decay  rates  and  diffusion  constemts  have  been  reproduced  exactly.  We  have  given  a 
very  detailed  analysis  for  the  two  well  case,  and  have  obtained  rather  complete  results  for 
the  time  dependence  of  the  interwell  transition  probability. 

The  organization  of  the  paper  is  as  follows.  In  Section  2,  we  set  up  our  system 
Hamiltonian,  define  our  problem,  and  derive  an  exact  path  integral  expression  for  the 
transition  probability.  In  Section  3,  we  maJce  a  systematic  expansion  of  the  expression, 
and  derive  a  simple  yet  general  result,  useful  for  later  explorations.  Section  4  is  devoted 
to  the  two  well  case,  where  the  effects  of  temperature  and  disorder  are  studied  in  great 
detail.  In  section  5.  we  study  the  case  of  low  temperature  and  large  site  energy  disorder 
for  an  array  of  wells.  In  section  6,  we  consider  the  translationally  invariant  case  at  low 
and  high  temperatures.  Finally,  we  discuss  and  conclude  in  Section  7. 
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2.  General  formulation 


2.1  -The  system  and  the  problem 

We  begin  with  a  discussion  of  the  Hamiltoniein  governing  the  motion  of  our  system. 
Consider  first  the  case  of  an  external  particle  moving  in  or  on  the  surface  of  a  crystalline 
solid.  A  particular  example  in  our  mind  is  a  hydrogen  atom  (or  its  isotopes)  moving  on 
the  tungsten  surface[12, 13.26,27].  When  the  lattice  atoms  are  fixed  at  their  equilibrium 
positions,  the  particle  sees  a  periodic  array  of  potential  wells.  We  assume  that  the  tem¬ 
perature  is  sufficiently  low  such  that  only  the  lowest  level  of  each  potential  well  can  be 
occupied.  This  is  appropriately  described  by  a  tight  binding  model.  Let  [x  >  be  the 
Wannier  state  corresponding  to  a  well  centered  at  x,  and  the  particle  tunnels  from  one 
well  to  another  according  to 

ifp|x>=  e|x> -f  |x  +  1>,  (2.1) 

i 

where  e  is  the  site  energy  and  A  the  tunnelling  energy  splitting,  and  1  is  a  vector  connecting 
nearest  neighboring  sites. 

The  motion  of  the  lattice  atoms  is  described  by  a  phonon  field  whose  Hamiltonian  in 
the  absence  of  the  external  particle  is  given  by 

Hl  =  (2.2) 

k 

where  aj^  creates  a  phonon  of  wave  vector  k  (a  mode  branch  index  is  omitted  for  simpler 
notations).  When  the  lattice  atoms  are  displaced  from  their  equilibrium  positions,  the  site 
energies  of  the  particle  also  change.  To  first  order  in  these  displacement,  the  site  energy  e 
cheinges  by 

=  +  aL  J,  (2.3) 

where  N  is  the  total  number  of  lattice  atoms.  The  actual  form  of  g(k)  depends  on  the  type 
of  interaction  between  the  particle  and  the  lattice  atoms.  [28]  In  an  ionic  crystal,  the  particle 
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interacts  strongly  with  the  optical  longitudinal  phonons  through  the  dipole  field  that  they 
generate.  The  coupling  function  takes  the  form  g{k)  -  n/k,  where  k  is  independent  of 
the  particle  mass.  In  a  piezoelectric  material,  the  particle  also  interacts  strongly  with 
the  acoustical  longitudinal  phonons.  The  coupling  function  is  y(k)  =  «:(k)/y^,  and  it 
is  again  independent  of  the  particle  mass.  In  Appendix  A,  we  give  a  derivation  of  (2.3) 
for  a  short  ranged  repulsive  interaction,  where  we  find  that  the  particle  interacts  with  the 
longitudinal  phonons  with  a  coupling  function  ^(k)  ~  k/ yjuk  for  small  k.  The  overall 
energy  scale  of  the  coupling  is  of  the  order  of 

In  this  expression.  St  is  the  energy  level  spacing  in  a  well,  M  the  mass  of  the  lattice  atoms. 

the  Debye  frequency  of  the  lattice,  and  I  the  lattice  constant.  The  coupling  strongly 
depends  on  the  particle  mass  through  the  level  spacing  St. 

If  (2.3)  is  the  only  form  of  particle-phonon  interaction,  our  system  Hamiltonian  is 
then  given  by 

Htot  =  Hp  +  Hl  +  Hi.  (2.5) 

Before  we  proceed  further,  we  would  like  to  remark  on  the  reality  of  this  Hamiltonian.  The 
reduction  of  the  original  continuous  system  to  the  tight  binding  form  has  been  worked  out 
in  various  ways.  Sethna[29]  has  given  a  quite  detailed  analysis,  showing  that  the  different 
reduction  schemes  are  equally  valid,  but  the  interpretation  of  A  may  be  different.  Our 
interpretation  of  A  as  the  bare  tunneling  amplitude  relies  on  the  assumption  that  the  intra 
well  level  spacing  is  large  compared  with  kpT  and  with  the  Debye  energy.  Second,  the 
hopping  amplitude  A  may  also  have  an  explicit  dependence  on  the  the  positions  of  the 
lattice  atoms  emd  therefore  giving  rise  to  tin  off-diagonaJ  coupling  to  the  phonons.  This 
is  the  so  called  ‘effect  of  fluctuation  preparation  of  barrier'. [30]  It  is  estimated,  however, 
that  this  effect  is  not  important  for  a  low  mass  particle,  such  as  a  hydrogen  atom,  moving 
on  a  heavy  substrate,  such  eis  the  tungsten  metal.  [31] 
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The  Hamiltoniaxi  (2.5)  is  translationally  invariant.  In  our  later  discussions,  we  will, 
however,  generalize  our  system  to  include  situations  without  translational  invariance:  1 ) 
the  site  energies  may  have  a  positional  dependence,  and  2)  only  a  subset  of  the  sites  may 
be  accessible  to  the  particle.  We  thus  modify  (2.1)  to 

^fplx>=  ex|x> -t-A  ^  |x  +  1>,  (2.6) 

I 

where  {x}  is  the  tight  binding  lattice  or  a  subset  of  it.  Due  to  the  smallness  of  the  tunneling 
energy  splitting,  static  disorder  due  to  impurities  or  strain  can  play  an  important  role.  .A.lso. 
the  inclusion  of  a  two  site  case  allows  for  a  very  detailed  analysis  of  the  fundamental  step 
of  the  particle  motion,  and  allows  for  direct  comparison  with  results  in  the  literature. 

We  wish  to  study  the  following  problem.  Suppose  the  particle  was  at  Xo  at  time 
t=0.  tmd  the  phonons  were  in  equilibrium  at  temperature  T  with  the  particle  fixed  at  xq. 
we  ask  for  the  probability  of  finding  the  particle  at  X/  at  later  times,  regardless  of  the 
final  states  of  the  phonon  field.  It  is  however  more  convenient  to  pose  the  problem  in  an 
alternative  but  equivalent  way. [19]  We  take  a  time  to  in  the  distant  past.  We  set  up  the 
initial  distribution  of  the  system  as 


P«o((<o)  =  |xo  Xxole  ^^‘/tr(e  (2.7) 

where  3  =  l/(kgT).  We  then  allow  the  system  to  evolve  according  to 

Htotit)  =  Hp{t)  +  HL+H[.  (2.S) 

where  Hp(t)  is  the  same  as  given  in  (2.6)  but  with  A  replaced  by  A^(t),  with  d{t)  being 
the  unit  step  function.  We  finally  take  to  — ►  -oo,  and  ask  for  the  probability  P(x/.Xo.t) 
of  finding  the  particle  at  Xf  at  time  t  >  0.  The  evolution  of  the  system  in  the  infinite 
time  interval  (to  <  t  <  0)  should  have  prepare  the  system  at  t  =  0  in  equilibrium  with  the 
particle  fixed  at  Xq. 
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2.2  -Path  integral  expression  for  P(x/,Xo,f) 

In  this  subsection,  we  will  derive  a  path  integral  expression  for  the  probability 
P(x/.Xo,t).  The  mathematical  details  will  be  presented  in  Appendices  B  and  C  and 
references  therein.  Here  we  just  illustrate  the  basic  ideas  and  establish  notations  used 
throughout  this  paper.  It  is  sufficient  to  know  the  reduced  density  matrix  p(t)  for  the 
particle,  defined  as  the  total  density  matrix  of  the  system  with  the  phonon  degrees  of 
freedom  traced  out.  In  the  absence  of  the  pairticle  phonon  interactions,  the  reduced  density 
matrix  satisfies  the  Liouvillian  equation 

^p{x.y,t)=  Cp{x.y;x'y'\t)p{x\y'j).  (2.9) 

x'.y' 

where  £p  is  the  Liouvillian  operator  corresponding  to  the  particle  Hamiltonian  Hp{t). 
defined  as 

£p(x.y;x',y';t)  =  ^{<x\Hp{t)\x' ><y'\y>  -  < x|x' >< y'lPp(t)|y  >},  (2.10) 

in, 

which  is  nothing  but  the  commutator  [Hp(t),  ■  \/iih). 

The  Liouvillian  equation  (2-9)  has  the  same  linear  structure  as  the  Schrodinger  equa¬ 
tion.  The  real  space  x  is  replaced  with  the  phase  space  (x,  y).[22]  The  wave  function 
is  replaced  with  the  density  matrix.  The  Hamiltonian  is  replaced  with  the  Liouvillian. 
Therefore,  ail  the  algebraic  techniques  developed  for  solving  Schrodinger  equations  can  be 
transplanted  here.  In  order  to  simplify  notations,  we  will  now  on  denote  the  pair  (x.  y) 
by  c,  so  that  p(z.t)  =  p(x.y.f)  and  Cp{z,z'J)  =  £p(x,  y;  x'.  y';  t ).  Also,  it  will  be  con¬ 
venient  to  introduce  a  ‘bra-ket’  notation  (not  to  be  confused  with  those  for  the  states  as 
in  (2.10))  for  the  quantities  in  the  linear  equation  (2.9): 

\p(t)>  =  ^/?(z)|2> 

£p(<)  =  £p(2,2',t)  <r'l  (2-11) 

»  -/ 

<  z\z  >  =  Szz'  —  ^xx'*^yy'- 
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Equation  (2.9)  can  then  be  written  zis 


—  \p't)>=  Cp{t)\p{t)>  . 


(2.12) 


(2.13) 


This  is  integrated  with  the  initial  condition  |p(fo)>=  l.:o>-  -^o  =  (xq.Xo).  yielding 

lp(t)>=^  /  f  I  dti 

fi=0  ^0 

X  )...£ft(n  )lzo  >. 

where 

£;,(f)  =  (2.14) 


with  £,  (time  independent)  corresponding  (as  in  (2.10))  to  the  site  diagonal  part  oi  Hp(t) 
and  Ck{t)  corresponding  to  the  hopping  part  of  Hp(t).  The  probability  of  finding  the 
particle  at  x/  is  then  given  by 

<Zf\p{t)>=^  /  dtn  /  dtn-i.:  /  dti  y 

n=0''(o  J>o  dto 

X  <  Z/|£h(t„)|2„_i  ><  Zn-l|£ft(^n-l  )l-n-2  >  < -1  l£/i(  ^  1  )  I -0  >  (2.15) 

n 

X  exp{^  A(r_,)(fj+1  -  tj)}, 

j=0 

where  Zj  =  {Xj,yj )  with  y/  =  X/.  and  X{zj)  =  (cx,  -  Cy,  )/{ih). 

The  expression  (2.15)  is  in  fact  a  sum  over  all  possible  paths  Zr  in  the  phase  space 
starting  at  zq  and  ending  at  Zf.  A  particular  path  is  specified  by  the  number  of  jumps  n. 
the  jumping  instants  {tj).  and  the  'positions'  between  the  jumps.  More  specifically, 
a  path  Zr  of  n  jumps  is  a  sectionally  constant  function  of  time  r.  given  by 


Zr  =  Zj,  for  tj  <  T  <  tj+i.  J=0.1....n.  (2.16) 

9 

where  tn+i  =  t.  The  weighting  factor  is  a  product  of  jumping  amplitudes  and  the  phase 
factors  accumulated  between  the  jumps  due  to  the  site  energies. 
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In  the  presence  of  interactions  with  the  phonons,  the  expression  (2.15)  is  modified  by 
inserting  an  influence  functional  as  (see  Appendix  B) 


n=0''‘o  •''o  “'‘o  j 


X  <r/|Zl/i(#ri)|.:^n  — iX-^n  — l|'^/i(^n  — i)|2ti  — 2>  <’2ll-C/,(ti)|co>  (2.1() 

n 

X  exp{^  \{zj)[tj+i  -  f_,)}exp($[c]), 


where 


m--  f  dr  f 


|g(k)i- 

Nri^ 


(e-kx,  _  _  ,)e-.wx.  _  _  3)^-“^*]. 


(2.18) 


The  influence  functional  fully  takes  care  of  the  effects  of  the  phonons  on  the  particle 


motion. 


With  some  tedious  but  straightforward  manipulations  (see  Appendix  C),  the  expres¬ 
sion  (2.17)  can  be  rewritten  eis 

<Zf\p[t]>=^  /  dtn  /  dtn-l...  /  dti  } 

n-odo  Jo  Jo  j 

X  <[C^|T/i|Zn  —  — ...  <iri|T/i|iro^  (2.19) 

n 

X  exp{A(ro)fi  -I- ^  A(c;)(f;+1  -  fj)}  JJ  Fjy. 

>=i 

In  this  expression  £/,  is  the  renormalized  hopping  Liouvillian  operator,  given  by 

A  V—' 

<  z\Ch\z'  >  =  -^  -b  Ixy'ly  >  -  <x|x'  xy'jy  +  1>} 

‘  (2.20) 

A  =Aexp{  — -  J(k)coth(/5ftu;k/2)|l  —  e'***'  |^}, 


(2.20) 


with  li  being  any  one  lattice  vector  among  those  of  the  nearest  neighboring  ones  {!}.  and 

Fjj'  is  an  influence  factor  given  by 
F jji  =  F{tj  —  tji ,  Zj,  Zj—i,  Zji ,  Zj'  —  I ) 

=  exp{-  ^  J(k){coth(j3fiu,‘k/2)cosu;k(f;  -  ty){Uj  -  Uj-i){UJ,  -  L^-i)  (2.21) 

k 

-  isintJkff,  -  ty){Uj  -  -  F, :_,)}}. 


The  other  quantities  in  the  above  expressions  are  defined  eis 


i(k) 

Vj 


|g(k)P 

Nihu^f' 

_gikx,  _  gjkyy  ^ 
^g.kx,  ^  g.kyy 


And  it  is  understood  that  U-i  =  V-i  =  Un+i  =  0. 

It  is  seen  that  the  hopping  energy  A  gets  renormalized  by  the  phonons.  The  renor¬ 
malization  factor  is  just  the  usual  Frank-Condon  factor  in  small  polaron  theory.  [1]  In  the 
next  section,  we  will  have  a  closer  exaimination  of  the  physical  contents  in  the  expressions 
(2.19)-  (2.21),  and  will  present  a  useful  approximation  scheme  to  evaluate  them. 
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3.  The  ‘rare  hopping’  approximation 


•  '  -Diagrajn  expansion 

In  the  last  section,  we  have  derived  a  path  integral  expression  (2.19)  for  the  probability 
of  finding  the  particle  at  X/  at  time  <  >  0,  eissuming  the  particle  wets  at  Xq  at  t  =  0.  and 
that  the  phonons  were  in  equilibrium  with  the  peirticle  fixed  at  xq.  The  effects  of  the 
phonons  on  the  particle  motion  come  in  two  ways:  1)  to  renormalize  the  hopping  energy 
A  by  the  Frank-Condon  factor  (2.20),  and  2)  to  dress  up  the  particle  propagation  by 
the  influence  factors  (2.21)  for  each  pair  of  hopping  events.  An  influence  factor  Fjy 
depends  only  on  the  time  difference  {tj  —  ty )  between  th.e  two  hopping  events,  and  on  the 
■positions’  immediately  before  and  after  each  of  the  hops.  From  now  on.  we  will  use  the 
expression  (2.19)  as  the  formal  definition  of  a  propagator,  G(z/,zo,t),  even  if  zq  and  z/ 
are  not  in  diagonal  form.  For  diagonal  {zf,zo),  G  coincides  with  the  propagator  of  the 
reduced  density  matrix. 

If  the  influence  factors  were  unity,  the  expression  (2.19)  would  be  equal  to 

<  Zf\e^^-+^''^*\zo  >=  Goizf,  -0,0,  (3.1) 

that  is,  the  particle  would  propagate  under  the  action  of  the  renormalized  particle  Liou- 
villian 

Cp  =  £,  +  £/,.  (3.2) 

The  typical  time  scale  between  two  hopping  events  is  ^/A. 

On  the  other  hand,  an  influence  factor  differs  appreciably  from  unity  only  in  a  typical 
time  scale  of  l/uj,  assuming  that  the  coupling  function  J(k)  varies  smoothly.  There  is.  in 
general,  a  long  time  tail  (decaying  to  zero  as  an  inverse  power  of  time)  in  F,y-l  associated 
with  the  low  frequency  end  of  the  phonon  spectrum,  but  this  will  be  small  due  to  the  small 
spectral  weight  there. 


13 


We  will  now  on  treat  i\/{hud)  as  a  small  parameter.  It  is  then  useful  to  write 


Fjji  —  1  +  {F jj'  1),  (3.3) 

and  to  expand  the  product  ]][  Fjy  in  (2.19)  in  powers  of  (Fjy  —  1).  This  is  analogous  to  the 
technique  used  in  the  low  density  expansion  of  thermodynamic  quantities  for  a  hard  core 
interacting  gas.  There,  the  expansion  is  in  terms  of  —  1  instead  of  the  potential  V'. 

Here,  we  may  say  that  we  are  dealing  with  a  ‘low  density’  (in  time)  of  scattering  events. 

Regarding  the  exp’-ession  (2.19)  as  a  propagator  G(z f,  zo,t),  then  it  is  equad  to 
Go{zf,  zq,  t)  dressed  with  zero,  one,  two,  ...factors  of  the  retarded  self-interactions  {Fjj-  —  1 ) 
in  all  passible  ways.  Diagramatically  this  is  shown  in  Figure  1,  where  a  heavy  directed  line 
represents  the  dressed  propagator  G,  a  light  directed  line  the  ‘free’  propagator  Go ,  and  the 
dashed  lines  the  self-interactions  (Fn'  —  1).  Two  self-interactions  can  share  at  most  one 
hopping  time.  The  detailed  structures  of  these  diagrams  will  be  made  clear  below. 
Consider  first  the  term  with  a  single  self-interaction, 


~  ft  ftn  fti 

Gl{Zf,Zo,t)=\]  /  dtn  /  dtn-l...  dti  } 
„_0-'0  Jo  Jo 


X  <  2^ jZji  — 1  — 1  l^/i  l^n  — 2  ^  -1  l^/i  I'^O  ^ 


(3.4) 


X  exp{A(2o)#l  +  ^  ^  (Fjy  -  1). 

j=i  ">;>;'>! 

The  structure  of  this  expression  is  more  transparent,  if  we  make  a  Laplace  transform  of  it. 


The  result  is 
Gi(zf,zo,  A) 


X{ 


oo  n  J  —  l 

-<=<=>  t.  =  2j  =  2j'  =  1{zj} 


z j\C,^\Zfi—\  ^ 


•...  <  2j'  +  i  |.C>i  jCj'  > 


A- A(2„-i) 


1 


1 


A-A(2/) 

"  - A-A(.,.)+.w. 

X  {"7  r"7  ■  ^  j' —  h\^  ]' —2  ^  T  w  T'"  |^/i  l^'O  ^  ^ 

A-A(2j-_i)  '  'J  \-\{zy-2)  A-A(2o) 

X  {<Zj\Ch\Z]-\  >  F{uJ,Zj.Zj-uZj>,Zj'-i)  <Zy\Ck\Zj'-i  >}, 


} 


(3.5) 
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where  F  is  defined  by 


/  (Lje  F(ui,Zj,Zj-i,Zj> 

J-oo  (3.6) 

=  F{  T,  Zj ,  Zj  —  I  ,  Zji ,  Zji  ^i)  —  1 . 

with  F(...)  given  in  (2.21).  we  now  change  the  labels  such  that  z*  =  Zj,  z^  =  Zj-\,  z^  —  -y  • 
and  2^  =  Zj'-i-  After  summing  up  the  intermediate  positions  other  than  these  c's.  the 
expression  (3.5)  becomes 


oo  n  7  —  1 


Gi{Zf,  Z 


n=2  7=2  7'  =  1  {i*} 


X  { <2*14 1--’ >  <2'*|4|--‘>). 


Now,  the  powers  n— j,j  —  j'  —  I  and  j'  —  1  are  remged  from  0  to  oc  independent  of  one 
another,  and  the  summation  over  them  gives 


Gi(2/,2o.  A)  =  f  dw  ^ 
J  —  OO  i  1 


1 


X  <-/It - ■? - - r - r  'I  ■  I 

X  <z*\Ch\z^>  F{^,z\z\z'^.z^)  <z^\Ch\z^>  . 


\  —  Cs  —  Ch  -‘r  i 


\-Cs-Ch 


Notice  that  the  quantity 


appeeu-ed  in  (3.8)  is  just  the  Laplace  transform  of  Go(^,^'.t)  defined  in  (3.1). 

As  a  result  of  the  above  analysis,  the  first  order  diagram  in  Figure  1,  tifter  Laplace 
transformation,  has  the  detailed  structure  shown  in  Figure  2.  A  directed  line  still  represents 
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a  free  propagator  Go,  with  its  position  variables  labelled  at  the  ends  of  the  line,  and  its 
Laplace  \’ariable  labelled  beside  the  line.  The  hatched  circle  with  four  legs  represents 
the  self-interaction  F,  with  its  position  variables  labelled  at  the  ends  of  the  legs,  and  its 
frequency  labelled  beside  the  circle.  Finally,  the  undirected  short  lines  between  pairs  of 
dots  are  the  hops,  <  >  and  <z‘^\Ch\z^  >■  The  ‘value’  Gi(z/,::o,  A)  of  the  diagram 

is  equal  to  the  product  of  the  quantities  associated  with  the  p^lrts  of  the  diagram,  with  the 
internal  positions  summed  up  zmd  the  frequency  integrated  out. 

Similarly,  the  second  order  diagrams  in  Figure  1  can  be  shown  to  have  the  detailed 
structures  displayed  in  Figure  3.  The  meanings  of  the  various  parts  of  these  diagreims 
are  the  same  as  for  the  first  order  one.  Once  the  first  and  second  order  diagrams  are 
understood,  it  is  not  difficult  to  imagine  how  the  higher  order  ones  should  look  like.  In 
a  mth  order  diagram,  the  otherwise  free  propagation  of  the  particle  is  interrupted  by  m 
self- interactions  between  m  pairs  of  hopping  steps.  Two  self-interactions  can  share  at  most 
one  hopping  step,  but  one  hopping  step  can  be  shared  by  any  number  of  self-interactions. 
Each  allowed  topologically  distinct  diagrtim  contribute  a  term  for  G. 

The  propagator  G  satisfies  the  Dyson  equation  (in  operator  notation). 

G(A)  =  Go(A)E(A)G(A),  (3.10) 

where  the  self-energy 

<r|E(A)|F  A)  (3.11) 

is  given  by  the  diagrams  shown  in  Figure  4.  These  diagrams  start  and  end  with  hops  and 
each  of  them  cannot  be  disconnected  by  removing  one  Go  line.  The  first  order  diagram 
has  the  value 

<z|Ej(A)|F> 

dw  F(u;,^,^l,^^,  r')  (3  12) 

X  <  z\Ch\z\  >  Gq{z\ ,  r'l ,  a  -h  tu;)  <  >, 

as  can  edso  be  read  off  from  (3.8). 
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3.2  -The  approximation 


We  then  maice  the  approximation 


(3.13) 


The  ignored  higher  order  terms  contain  at  least  three  ‘hops'  and  two  self-interactions, 
and  therefore  they  axe  at  least  A/ihu/d)  times  smaller  than  the  kept  one.  With  this 
approximation,  we  can  write 


G(A)  = 


1 

A-£p-Ei(A)- 


(3.14) 


The  Laplace  transform  of  the  probability  P(x/,Xo,<)  is  then  given  by 

P(x/,Xo,  A)  =<2/1G(A)1xo  > 


=  <2/| 


A-£p-!:i(A) 


Fo  >• 


(3.15) 


where  z/  and  cq  are  in  diagonal  form.  Thus  the  problem  becomes  solving  for  Ili(A),  finding 
the  inverse  ofA  —  £p  —  Ti(A)  and  inverting  the  Laplace  transform. 

The  expressions  (3.12)  and  (3.15)  are  the  central  result  of  this  section.  The  operator 
inversion  is  facilitated  by  studying  the  eigenvalue  problem. 

A|i/;>=  [£p  -I-  Ti(A)]  ii:;>  . 


This  is  amalogous  to  the  generalized  Boltzmann  equations  in  the  literature.  [23-25]  The 
problem  is  further  simplified  by  restricting  to  small  values  of  A.  if  we  are  interested  only 
in  the  long  time  behavior  of  the  particle. 

The  rest  of  this  subsection  deals  with  the  detailization  of  the  e.xpression  (3.12).  The 
reader  may  wish  to  skip  this  part,  and  come  back  later  when  the  formulas  are  referred  to. 

It  will  be  useful  to  write  (3.12)  in  a  more  detailed  form.  It  is  easy  to  see  that  (3.12) 
is  the  Laplace  transform  of  the  following  expression: 


<r|!:i(r)(r'>=  [F(r.  r,  n c')  -  1]  < X  d  >< 1 A I-' >  •  (3-16) 
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From  (2.20)  we  have. 


^  {^x,xi  +  l^y,yi  ^  y  I  ^x.xi  <^y,yi— 1  <  X.  y  +  l|  }. 

1 

l^lX^jl^^l^  >=  ^  ^  "  {^x' ,x'+t'^y',  ,y' 1^  "t"  ^  '  y  ^  ""^Xj  ,x' ‘^y'j  ,y' —  1' |x  .  y  —1  >}. 

r 

(3.17) 

Substituting  this  into  (3.16)  and  using  the  definition  (2.21)  for  F{t,z,zi,z\,z').  we  have 
<z\^,{t)\z>> 

=  {j?  5]{-  <X  -  Lyle^^Hx'  +  l',y'  >  Fn-(r  +  ihdl2,^-  x'  -  1/2  -  172) 

i.r 

-  <x,y +  lle^'’nx',y' -l'>  Fii<(r-z;id/2,y-y'  +  l/2  +  l72)  (3-18) 

+  <x  -  l.y|e^'*7x',y'  -  l'>  Fn.{T  -  iht3/2,x-  y'  -  1/2  +  172) 


+  <x.y +  l|e'^'’7x'  + I', y'>  Fu-(r  +  ead/2,y-x'  +1/2-172)}. 


where 


Fndr.x)  =  exp  ^-4^ 


J(k)sin(kI/2)sin(kr/2)  .^x 

— si;ih(d^72) — 


(3.19) 


The  Laplace  transform  of  (3.18)  then  gives 


<rlTi(A)lr'> 

=  <x-lyl 


i.r  >*/.k 


g;JAu//2g.k(x-x'-l/2-l72) 
\  +  lU!  —  Cp 


x'  +  l',y'> 


-  <x,y  + 1|- 


+  <x  -  l,y|- 


+  <x.y  +  H 


^-iiftw/2„.k(y-y'+l/2+l72) 

^ ^ ? - Ix'.y'  -1': 

\  +  lU-Cp  ' 

„  —  jcJftw/2pik(x-y' -1/2  +  172) 

^ - |x', y'  -  I'  > 

X  +  iu!  —  Cp 

^  +  J/lw/2„ik(y-x'  +  l/2-l72) 

^  ^  \..i  ,  1/ 


x.y'-l'>  (3.20) 


A  +  luj  —  Cl 


lx'  +  l',y'>}. 


where  Fud^.k)  is  the  Fourier  transform  of  Fn'('''«x)  given  by 
5]Fudu;.k)e-‘""e'‘‘’‘ 


=  exp 


_  J(k)sin(kl/2)5in(kl72)  ^ 

+  sinh(/Jftuk/2)  '  '  / 


(3.21) 
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At  temperatures  low  compared  with  the  Debye  temperature,  the  exponent  in  the  above 
exponential  is  small.  To  first  order  in  the  exponent,  we  have 


/(k)  sin(kI/2)  sin(kIV2) 
sinh(/3^t4;k/2) 


(3.22) 


It  is  not  difficult  to  see  that  we  can  obtmn  the  same  result  from  (3.16)  if  we  replace 
F{t.z,Zx.z[,z')  -  1  by 

-  J (k) {coth( f3fiui,/2)  cos uuiT)iU  -  Ux){U["  -  U'*) 

k 

~  is\nuu{T)[U  -V* )), 

which  is  the  argument  of  the  exponential  in  the  definition  (2.21). 


:3.23) 
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4.  Two  level  system 


In  this  section,  we  consider  the  case  that  the  particle  is  restricted  to  only  two  sites, 
which  axe  located,  for  instance,  at  x  =  li/2  and  x  =  — li/2.  Our  system  then  becomes  a 
two  level  system  coupled  to  a  heat  bath  of  phonons.  The  analysis  will  be  done  in  three 
subsection:  1)  general  considerations,  2)  the  case  of  low  temperature  but  arbitrary  bias, 
and  3)  the  case  of  zero  bias  but  arbitrary  temperature. 

4.1  -General  considerations 
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The  spectral  density  J(uj),  with  a  dimension  of  frequency,  contains  all  the  information 
of  the  phonons  relevant  to  the  dynamics  of  the  two  level  system.  [19]  It  vanishes  beyond 
a  frequency  around  a,(f.  At  the  low  frequency  side,  it  vanishes  as  a  power  of  For  a 
short  ranged  repulsive  interaction  between  the  particle  and  the  lattice  atoms.  l3(k)|^  goes 
to  zero  linearly  with  .*;w.  so  J{u!)  vanishes  like  at  low  frequency.  For  a  piezoelectric 
coupling  Jiui)  vanishes  like  In  Ref.  19,  the  genereJ  case  of  J(u;)  ~  has  been  studied 
for  the  low  temperature  dynamics  of  a  two  level  system.  Our  case  lies  well  in  the  regime  of 
superohmic  coupling  (s  >  1),  where  the  low  spectral  density  at  small  uj  makes  the  coupling 
effectively  weak  at  low  temperatures.  However,  the  low  frequency  behavior  of  J{uj)  will 
not  be  so  important  for  higher  temperatures. 

In  order  to  find  the  probability  P(x/,Xo,A)  (3.15),  we  need  to  know  Ei(A)  which  is 
the  Laplace  transform  of  Ei(0  given  by 

<z\Z,{t)\z'> 

=  '^{F(t,z,zuz[,z')  -  1)  (4  6) 

X  <r|£hl2l  ><2i|e^'’'|2'l  ><2'J£ftl2'>  . 


Throughout  this  section  we  will  use  z  to  denote  the  pair  (x,y).  Then  the  operators  £/, 
and  Cp  are  given  by 

Cp  =  £/,+£, 

A 

£/,lx.y>  =  — (|x.y>  -|x.y>)  (4.7) 

in 

£,|x,y>  =  -  y)|x.y>. 

in 

where  x  =  — x.  y  =  -y,  and  ex  is  the  site  energj’  of  the  site  x  =  xli/2.  Therefore, 


<.r,  y|<5rii<5yiy  ) 


Iz'i  ><2'i|£^|2'>  =  — (|x',y'>  -  [x'.y  >  dp^pSy.^g,). 


(4.8) 
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Substituting  the  above  relations  into  (4.5),  we  have 

h 

-  <x,y\e^^^\x\y'>  -  1) 

+  <x.y|e^^'|x'.y'> 

+  <x,y|e^'’‘|x',y'>  _  i)}. 

where 

du 

git  ±  ih3/2)  =  4  I  ~J{uj){coth.{3hui/2)cos{u;t)^^isin(u.'i)] 

Jo 

du  J{^) 

Jq  u;2  sinh(/?rtu;/2) 

We  now  introduce  the  Pauli  matrices  {(7^,/i  =  1,2,3},  defined  as 

(Ti|x,y>=  |x,y> 

(y2\x,y>=  ix\x,y> 

(T3lx.y>=  x|x,y>  . 

The  particle  Liouvillian  Cp  can  then  be  written  eis 


(4.10) 


(4.11) 


Cp  =  ^[(^CTi  +  ecTa)  -  (Acr'i  +  ecr^)], 


(4.12) 


where  a'^^  acts  on  the  y  variable  the  same  way  as  cr^  on  the  x  variable.  It  can  then  be 


shown  that 


where 


::i(i) 

=  ( ~  {-(T\e^’’^aiC{t  +  ih3/2)  -  cr3crie^'’'(7icr3  5(t  +  ih3/2) 
h 


-  a[e^^'cr[C{t  -  ih3/2)  -  (7>; cr;5(f  -  ih3l2) 

+  -  ih3l2)  -  o-3<7ie^'’'<T;<7;S(t  -  «7id/2) 

+  G\e^‘‘*axC{t  +  ih3l2)  -  a'^(T\e^^^a^(T:iS{t  +  ih3l2)], 

C(t)  =  cosh(y(t))  -  1 
Sit)  =  sinh(y(0)- 


(4.13) 


(4.14) 
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The  operator  £p  is  diagonalized  by  the  following  transformation 


in 


(4.15) 


where 


n  = 

A 


sin(20)  = 


a 


(4.16) 


cos(29)  =  -. 

Under  the  same  transformation,  the  self-energy  operator  Si  (A)  becomes 


i(A) 


=  — (  — )^y  du,'[C(w)Ai  -h  5(w)A2], 


(4.17 


where  C{uj)  and  S{u>)  are  the  Fourier  transform  of  C{t)  and  S{t)  respectively,  and  -4i  and 
,42  are  given  by 


e  A 

(-(7i  +  -^3) 


/  ^  I  \ _ 

x  +  iu  +  i{<73 '  n 


I  ^  ^  I  \ 

+  77<^3) 


^  —  Jhw/2 


£  A  _ 


I  ^  ^  (  N 


.  €  ,  A  , 


„^hw/2 


Q''^'  \  +  i^'  +  i{a^  ^ 


.42  =  {<73<^1 


^JhiA//2 


.  f  f 

+  cr^a 


A  -)-  ^u;  -b  )((T3  —  (7'^)Q/h 

g-Jftw/2 

^  A  "b  lU!  F  i( ^3  —  *^3  ^ 


(71(73 


-b  (T3<Ti 


.  /  f 

-b  cr.cr. 


g-JAu//2 


A  -b  tw  +  i((T3  —  a’^)?tlh 

g3ftw/2 


I  I 
(7i(73 


(7i(73 


<yicr3}- 


(4.1S) 


^  ^  X  +  iu  +  iicr^  —  (T'^)Clfh 
In  terms  of  the  transformed  operators,  the  probability  (3.15)  can  be  written  as 

1 


P(x/,xo,A)  =<  X/. 


A  -b  j(cT3  -  (7^)n/^  -  Si(A) 


e'®('^’+'^^)lxo,xo>  .  (4.19) 
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When  X/  =  — l.xo  =  +1.  we  have 


P(-,+.A) 


A  +  i((73  —  a'.^)Q/h  —  Si{A) 


1^’++  >' 


14.20) 


where 


K’++  >—  .^[(1  +  ^l++  ^  +(1  —  Q  )|- 


A,  A, 

->+||+->+||_+>| 
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4.2  -Low  temperature  and  arbitrary  bicis 

When  the  temperature  is  low  compared  with  ♦^'^e  Debye  temperature,  the  function 
g{t)  (real  time)  in  (4.10)  will  be  small  di.e  to  the  high  frequency  cut  off  by  sinh(  J^u,V2) 
and  the  low  spectral  weight  of  j{uj)  at  low  frequencies.  We  cem  therefore  set  C{t)  =  0  and 
S{t)  =  g{t)  to  leading  order  in  ksT/ihud)-  The  expression  (4.17)  then  becomes 

2J{u;) 


A  ,  du; 
Ti(A)  =  -(-)2y 


(4.22) 


h'  sinh(3^w/2)' 

where  J{u;}  =  — J(— u.')  for  <  0. 

Now.  both  3]i(A)  and  i{crj  —cr'-^)Q/h  are  block  diagonal  in  the  sectors  with  <73(73  =  ±1. 
Written  in  matrix  form.  A  4-  <(<73  —  a'j)Q/fi  —  Di{A)  is  equal  to 


A  4-  a( A)  —  A) 
— a(  A)  A  +  6(  A ) 


in  the  subspace  of  {I+4.  >.  |__  >}.  and  it  is  equal  to 


\ -i- i2Q/h  +  c(\)  c(A) 

c(A)  X  —  i2Q/h  +  c(X) 


(4.23) 


(4.24) 


in  the  subspace  of  — (->),  where 


a(X) 

A  . 

du! 

L 

u}^ 

A  , 

,00 

du! 

6(A) 

=  ] 

.2 

n  J 

1 

—  00 

tjj 

c(A) 

A  , 

fCO 

du! 

6  j 

L 

u!^ 

2J(c 


2Ae-‘^ftw/2 


2^(-) 


2Ae4ftu(/2 


(4.25) 


2J[^)  2Xe-^^^^^ 
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The  matrices  (4.23)  and  (4.24)  can  be  easily  inverted  to  give  an  explicit  expression 
for  P(  — ,-)-,A)  in  (4.20)  as 

_ 


P(_  ,  A)  =  -  -  n(^(^)-^(A))+(^)^A  _  1  _ 

2A  2A(A +  a(A) +  &(A))  2  A2  +  2Ac(A)  +  4f)2/^^  ' 

This  expression  has  four  poles  located  at 

A  =  0. 

A  =  -(a(A)  +  6(A))  ^  -(a(0+)  +  6(0+)), 

A  i2Q/n-c(0+  +i2Ctlh), 

A  25  —i2Cl/h  —  c(0+  —  i2Q,/h), 

respectively.  With  their  residues  approximated  to  first  order  in  (A/ 6)^.  we  have 

^(-.+,A)=i[l  +  ^tanh(/9fi)] 

^tanh(;3Q)  +  (^)^ 

2(A  +  2r) 

A  +  r 


(4.26) 


(4.2'; 


(4.2S) 


where 


-( —  f - 

2^fi^  (A  +  r)2  +(2fi'/6)2’ 


r  =  7r(^)^/(2fi/6)coth(dn) 

4  c 

O'  -  O  a- ^^^2  8J(u;)coth(i6u;/2)^ 

\  I  u;2  (2fi/6)2-+2  ) 


(4.29) 


Recall  that  Q.  is  the  eigenvalue  of  the  renormalized  particle  Hamiltonian.  Hp.  now  Q  is 
further  renormalized  to  fl'  by  the  residual  interaction  with  the  phonons.  Afte’-  transforming 
back  to  the  time  domain,  we  finallv  have 


P(-,  +J)  =i[l  +  ^  tanh(.3n)] 


-  ^(^tanh(^Q)  +  (^)^]e 

cos(2n'f/6)e-^‘. 


(4.30) 


This  agrees  with  the  results  of  Leggett  et  al[19]  in  the  limits  of  zero  bizis  (e=0)  and  large 
bias  (e  >>  A),  and  it  fills  out  the  intermediate  region  where  the  non-interacting  blips 
approximation  used  by  these  authors  is  invalid. 
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The  above  result  may  be  interpreted  as  follows.  At  low  enough  temperatures,  our 
system  is  effectively  described  by  the  renormalized  particle  Hamiltonian  Hp  =  eay  +  A<7i. 
which  is  weahly  perturbed  by  the  residual  interaction  with  the  phonons.  [7]  In  the  repre¬ 
sentation  where  Hp  is  diagonalized,  the  diagonal  parts  of  the  particle  density  matrix  rel5LX 
exponentially  at  a  rate  of  2r  to  the  equilibrium  distribution  exp{—i3Hp)/tT  exp{  —  3Hp), 
while  the  off-diagonai  parts  oscillates  with  a  frequency  20.' /h  %  20/h  and  damps  to  zero 
at  a  rate  of  F,  the  inverse  life  time  of  the  eigenstates  of  Hp.  The  probability  P(  —  .  -l-.t) 
for  the  particle  to  go  from  site  x  =  -|-1  to  site  x  =  —  1  ,  will  have  contributions  from  both 
the  diagonal  and  off-diagonal  parts  of  the  density  matrix.  In  the  limit  of  e  >>  A.  the  new 
representation  is  almost  site  diagonal,  so  that  P(  — .  +,t)  mainly  comes  from  the  diagonal 
parts  and  relaxes  purely  exponentiailly.  In  the  other  limit,  e  <<  A,  we  should  have 

~  —  p(— '  +>  0] 

where  p  is  the  density  matrix  p  in  the  new  representation,  and  we  have  used  the  fact  that 
the  trace  of  a  density  matrix  is  unity. 

4.3  -Zero  bia.s  and  arbitrary  temperature 


When  the  potential  bias  e  is  zero,  A  -I-  -  (T''i)0/h  -  !:i(A)  is  again  block  diagonal 

in  the  sectors  with  <73cr3  =  ±1.  In  matrix  form,  it  is  equal  to 


A  -h  a(  A)  — i(A) 

— a( A)  A  -I-  6(A) 


(4.31) 


in  the  subspace  of  and  it  is  equal  to 


A+  -i-  d(  A+ )  -f-  c(  A )  c(  A ) 

c(A)  A_-|-<i(A_)-(-c(A) 


26 


(4.32) 


in  the  subspace  of  {!  +  -  >,  1 >},  where  A±  =  A  ±  i2A/h.  eind 

d(X)  =(^f  n 

^  ^-oo  A  +  za; 

With  the  above  matrices  inverted  and  substituted  into  (4.20),  we  have 


(4.33) 


A  +  M<i(A+ )  4-  d{A_ )] 


(4.34) 


where 


(4.36) 


2A  2  [A+  4-  (i(A4.)][A_  4-  d(A_)]  4-  2c(A)[A  4-  j(d(A+)  4-  d(X-  ))] 

Using  the  definitions  of  S(u;)  and  C(a’),  we  can  rewrite  c(A)  eind  d{X)  explicitly  as 

A 

c(A)=2(  — )^  /  sinh(/?(0)  cos(/(t)). 

A  (-i-So) 

d{X)  =  A{-f  /  (ite-^‘[cosh(i?(t))cos(/(t))  -  1], 

”  Jo 

where 

/■'*  da> 

i?(f)  =  4  /  — 7(u;)coth(d^u;/2)cos(u;t), 

Jo 

The  expression  (4.34)  is  valid  for  any  temperature.  In  the  last  subsection,  we  have 
considered  the  case  of  low  temperature,  where  we  found  an  underdamped  oscillation  for 
P(  — ,4-,t),  indicating  a  coherent  motion  of  the  particle.  Here  we  explore  the  possibility  of 
a  transition  to  incoherent  motion  as  the  temperature  is  raised. 

We  first  give  an  estimate  of  the  poles  in  the  second  term  of  the  expression  (4.34).  So 
long  as  A  and  X±  axe  small  compared  with  aJj,  the  frequency  scale  of  R{t)  and  /(#),  we 
may  replace  them  by  zero  in  c(A)  and  d(X±).  We  then  find  that 


A  =  -[d(0)  +  c(0)]±\/[c(0)]2-4(-)2. 


(4.37) 


It  is  seen  that  the  poles  become  real  when  c(0)^  ^  4(  or  when 


4  r 

«  Jc 


dt  sinhiRit))  cos( lit))  >  1, 


(4.38) 


which  can  be  satisfied  only  when  the  temperature  is  high  enough.  The  transition  temper¬ 
ature  may  be  estimated  by  a  steepest  descent  evaluation  of  the  time  integral. 


dt  sinh( R[t))  cos(/(f )) 


if-.- 

\Jo  SI 


8  \7r  Jo  sinh{3fiu/2)J  \Jo 

The  left-hand  side  of  the  inequality  (4.38)  then  becomes, 


AJ(uj) 

sinh(/3?k<;/2) 


(4.39) 


14  (I  r 

8  h  \7r  Jo 


1  A  /2fcBT 

8  ^  V 


sinh(3^/2) 

r  —Ji^) 

Jo 


-[  sinh(/3 


(4.40) 


where  in  the  last  step  we  have  made  a  high  temperature  expansion.  This  result  can  also  be 
obtained  from  the  condition  that  the  life  time  of  bonding  or  anti-bonding  state  be  shorter 
than  2^,  a  time  scale  associated  with  the  level  spacing  2 A. 

Right  at  the  transition  temperature,  the  two  poles  merge  together  at 

\  =  -[d(0)  +  c{0)] 

-  ^  f^^-To/T 

^  4V2ftfcBrVro 


(4.41) 


where 


ro  =  7 —  /  — 

Jo 


J  ( u,’ ) . 


and  we  have  used  again  the  steepest  descent  method  and  high  temperature  expansion.  It 
is  seen  that  the  pole  is  of  order  A/ft  or  smal’er,  which  is  much  smaller  than  justifying 
the  replacements  of  A  and  \±  by  zero  in  c(A)  and  d(A±),  respectively. 

When  the  inequality  (4.38)  is  well  satisfied,  we  can  write  (4.34)  as 

1  1  ^+^(0) 


J_  _  1 _ A  +d(0) _ 

R(-,+.A)  _  2(A-hd(0)][A4- (i(0)  +  2ci0)] 

_J_  _  1 

~2A  2(A-f-r)' 


(4.43) 
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and  therefore 
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5.  Strong  site  energy  disorder  and  Markoving  motion 


In  this  section  we  consider  a  general  multi-site  case  with  site  energy  fluctuations  Se 
Icurge  compared  with  the  renormalized  hopping  energy  A.  For  the  two  site  case  studied 
in  the  last  section,  we  found  an  exponential  time  dependence  in  P(— suggesting 
a  Markoving  description  of  the  particle  motion.  This  motivates  us  to  derive  a  master 
equation  for  the  particle  motion  for  the  general  multi-site  case  here. 

We  consider  first  the  self-energy  operator  given  in  (3.12).  We  will  restrict  our  attention 
to  low  temperatures  (kgT  <<  hu’j)  to  simplify  calculations.  According  to  the  arguments 
given  at  the  end  of  Section  3,  we  may  make  the  approximation 


F(r,z,zi,2'i,z')  -  1 


J(k) 


smh(Jh(M'k/2 


{coswklr  -f  ih3/2)(e'^^  - 


+  coswk(r  -  ih3/2)(e‘'^y  -  ) 

-  coswk(r  -  1^3/2)16''^’^  - 

-coswk(r-t-tAd/2)(e‘*‘y -e-’*’'')  } 
Substituting  this  into  (3.16),  we  have  (in  operator  notation) 


+  cosu;k(r  -  e^J/2)[y.£p]e^'’l£p.y*] 
-cosu,'k(r  -  i7iJ/2)[X£p]e^'*^[£p.y*] 
-cosu;k(r  -h  jfid/2)[r.£p]e'^^^[£p,-Y*]  }, 


(0.1) 


(5.21 


where  [  .  ]  denotes  a  comutator,  and  A'  emd  K  are  operators  coreponding  to  and 
respectively. 

Suppose  |m  >  is  an  eigenstate  of  ffp  with  energy  €m,  then  an  eigenstate  of  £p  may 
be  written  as  jm.n  >  with  eigenvalue  (e^  —  Cn)/(i^)-  The  ‘diagonal’  states  {|n.n  >} 
?ire  special;  they  share  the  same  eigenvalue  0  of  Cp,  and  are  separated  from  others  by 
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a  frequency  scale  of  Stih.  Under  the  condition  St  »  A,  they  are  effectively  not  mixed 
with  the  others  by  !Ei(A).  The  inverse  in  (3.15)  may  therefore  be  carried  out  within  the 
subspace  of  the  diagoned  states. 

The  matrix  element  of  Di(r)  between  two  diagonal  states  is  given  by 
<  n,  n|Ei(r)ln'.  n'  > 

-7(k) 


=  E 


^  ft'' sinh(^ftu,'k/2) 

{coscJk(r  +  ih3l2)6nn'  ^  I  <nlXlni  >  l^(e„,  - 


ni 


:5.3) 


+  cosu.’k(r  -  ih3/2)6nn'  ^  1  <n2|y|n>  |^(e„  -  ""2) 

nj 

—  cosuJk(r  —  th3/2)  <nj,Y|n'><n'|y*|n>  (e„<  — 

-  cosu,’k(r  -h  ihd/2)  <  n'lF|n  ><  n|X"|n'  >  (€„  —  Cn' 
where,  with  the  definition  <x|n>=  ?Pn(x),  we  have 

<n\X\n' 

X 

<nlF|n'>=  ^Kiyy^^^n'iy). 
y 

We  then  Laplace  transform  (5.3)  and  set  A  =  0+.  yielding 
<  n.  njHif  A  =  0+)|n',  n'  > 

—  ^nn'  |<  n|Xlni  +(eni  —  en)/ft )  +  ~  (fni  —  f  n)/ft  ']}• 


(5.4) 


0.0 


The  density  matrix  then  satisfies  the  master  equation 


— p(n,n,f)  =  ^  '''^/^p(n',  n',  t)  -  *''' ^/^p(n.  n.  t )],  (5.6) 


where  p(n.n,f)  =<n,n\p{t)>,  and 


(5.7) 
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The  master  equation  involves  only  the  propabilities  on  the  eigenstates  oi  Hp.  indicating 
the  incoherence  of  the  particle  motion. 

The  physical  picture  is  very  clear.  When  the  site  energy  disorder  6e  exceeds  the 
renormalized  band  width  A.  the  particle  localizes  essentially  on  a  lattice  site.  The  particle 
is  mobalized  only  because  the  scattering  of  the  phonons,  which  knock  it  from  a  localized 
state  to  another.  This  picture  is  very  much  the  same  as  for  electron  hopping  conduction 
in  impurity  bands,  except  for  the  degeneracy  and  Fermi  statistics  in  the  latter  case. [32] 

The  master  equation  (5.6)  describes  a  random  walk  in  a  random  environment. [33] 
V’arious  techniques  have  been  developped  for  the  evaluation  of  the  diffusion  constant.  If 
the  hopping  probabilities  do  not  vary  too  much,  we  czm  ignore  the  percolation  effect. [34] 
The  diffusion  constant  may  be  estimated  by 

i 

i 

D  =  Wl\  [0.8]  1 


where  /  is  the  lattice  constant,  and  IT  the  typical  value  of  the  quantity  in  (5.7)  for  a  pair  of 
states  localized  at  nearest  neighboring  sites.  Here,  we  have  assumed  that  the  eigenenergy 
differences  lc„-  -  e„|  are  small  compared  with  kpT.  The  overlapping  factor  |  <  n|e‘‘‘'‘ln'  >  p 
can  be  czdculated  perturbatively  as  L  -1)1^-  where  1  is  the  vector  between  the 

sites  at  which  the  two  states  are  localized.  We  then  have 

”  sinh(JftWk/2) 


=  7t( 


h}  f-n  -  sinh(  J  lc„  -  en'l/2)  iV 

A  ,  T  knT 


^  -  1|2<5(u,-w  -  Icn  -  e.'l/^)  (5.9) 


-  jfn  -  €n'\/h). 


The  diffusion  constant  thus  depends  quadratically  on  the  renormalized  tunnuling  en¬ 
ergy  splitting  A,  linearly  on  temperature,  and  quadratically  on  the  particle  phonon  cou¬ 
pling  energy  ^(k).  The  dependence  on  the  strength  of  disorder  Se  differs  for  different  type 
of  couplings.  For  short  ranged  interaction  we  have  ]^(k))^  ~  (see  the  remarks  on 
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(2.3)),  therefore  D  ~  (^e)^.  The  diffusion  constant  increases  with  disorder!  The  more 
surprising  thing  occurs  for  a  piezoelectric  coupling,  for  which  |g(k)P  ~  — .  therefore 
is  independent  of  —  €„<.  The  diffusion  constant  is  independent  of  disorder!  Of  course, 
these  arguments  are  valid  only  in  the  regime  oi  A  <<  6e  <<  kgT  << 
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6.  Lattice  with  constant  site  energies 

6.1  -General  consideration 


In  the  case  of  translational  invariance,  it  will  be  convenient  to  work  in  the  momentum 
space.  Introducing  the  momentum  eigenstates, 

,-.pXe-.qyix,y>,  (6.1] 


X.y 


where  is  the  total  number  of  sites,  then  the  probability  (3.15)  may  be  written  as 
P(x/,Xo,  A) 

1  ,  ,  ,  (6.2) 


_  ^  g'((p'+q')xo-(p+q)x/] 


A-£p-Si(A) 


|p'.q'>  . 


p.q.p'.q' 

We  now  miike  a  variable  change  (p.  q)  — ►  (u,  v),  where  u  =  p  +  q  and  v  =  (p  —  q)/2.  then 
both  Cp  and  !Zi(A)  are  diagonal  in  u  (see  below),  and  we  may  write 
P(x/,xo.A)  =  ^^P(u.A)e*“<’‘'>-’‘^>. 


Ns 


1 


(6.3) 


V,V' 


A  -  £p(u)  -  i:i(u.  A) 


V  >  . 


In  general,  there  is  a  pole  A(u)  of  P(u,  A),  which  approaches  zero  as 


A(u)  =  —Du^  (6.4) 

for  small  u.  The  coefficient  D  is  the  diffusion  constant.  Therefore,  we  only  need  to  study 
P(u.  A)  for  small  u  and  A. 

We  first  derive  the  operator  Li(u,  A).  From  (3.20)  and  the  definition  (6.1).  we  have 


<p  q| 


i(A) 


=  < f  Z E “Ha < P- 


_g4ftw/2g-i(p-k/2)(l  +  l') 


i.r  w,k 


g-a)iu-/2g.(q-k/2)(l+I') 


<p.ql 


A  +  iu;  +  i(Cp  —  Sq— k)/^ 

.-A?ftu;/2--i(p-k/2)l+i(q+k/2)l' 

+  — .  ,  <p-k,q  +  k| 

A  -h  4-  z(ep_k  —  fq)/^ 

g4Au;/2g-i(p+k/2)r  +  i(q-k/2)l 


[6.5) 


-!-■ 


A  +  la;  +  j(€p  -  Cq-k)/^ 


<p-|-  k.q  -  k|  }. 
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where  Cp  is  the  band  energy  function  of  Ck-  It  is  seen  that  u  =  p  +  q  is  conserved  under 
i;i(A).  Within  the  sector  of  a  given  u,  we  write  <  p,ql  =<  v],  and  define  £p(u)  by 
<vl£p(u)  =<p,  q|£p.  therefore 


<  V|£p(u)  =<V|— (€v+u/2  -  ev-u/2)i 


(6.6) 


where  we  have  used  the  fact  that  Cp  =  c_p.  Similarly,  we  define  Si(u,  A)  by  <  v|Ei(  u.  A)  = 
<p,q|Ili{A).  We  therefore  have 


<  vj  Ei(u,  A) 


<v|  +  /2  <  V  +  k| 


(6.7) 


l,r  u;.k 

where  the  f  s  axe  defined  as 

g-.u(l  +  l')/2 


h  = 


/2  = 


g«ua+i')/2 


A  —  iu;  +  <(€v+k+u/2  ~  «v-u/2)/^  A  +  ZW  —  t(«v+k-u/2  “  ^v+u/2)/^’ 


,-iu(l-I')/2 


+ 


g.u(I-l')/2 


(6.8) 


A  —  zuz  +  i(ev+k+u/2  —  Cv-u/2)/^  A  +  zu;  -  i(ev+k-u/2  ~  fv+u/2)/^ 


and  we  have  used  the  fact  that  Fn<((*',k)  is  invariant  under  the  change  of  sign  in  k  and 
In  the  limit  of  u  =  0  and  A  =  0.  we  have 


<v|  !:i(0.0)  =  27r(^  k)e-'(v+k/2)(l+n^(^  -  (ev+k  -  fv)/^) 


i.r  u/,k 

X  <  v|  +  <  V  +  k|  ]. 


(6.9) 


6.2  -Low  temperature 

According  to  (6.3),  the  poles  of  P(u,  A)  are  given  by  the  eigenmodes  of  the  equation 

A|zi)>=  £p(u)lz/>> -f-Difu.  A)|z/j>  .  (6.10) 

We  are  interested  in  solutions  with  small  A  and  u.  We  can  set  Ei(u.  A)  =  Ei(O.O).  if  its  A 
derivative  is  small  compared  with  unity,  and  if  its  u  derivative  is  small  compared  with  the 
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typicad  band  velocity  (the  u  derivative  of  £p(u),  to  be  precise).  When  the  temperature  is 
not  high  compared  with  the  Debye  temperature,  these  conditions  are  satisfied.  Equation 
(6.10)  can  then  be  written  in  the  form  of  a  Boltzmann  equation 

dev 

X  lA’(v)  =u  — —  Tp(v) 
in 

v' 

where  Ww'  =  B^v'v  is  given  by 

Wvv'  =  27r(| )2  ^  ^  F„-(u;.  v'  -  v)e-<'''+'')''+'')/2^(u;  -  (6^-  -  ev)/^)-  (6.12) 

1,1' 

The  left  hand  side  of  (6.11)  corresponds  to  the  time  derivative  term.  The  first  term  on 
the  right  hand  side  of  the  equation  corresponds  to  a  drifting  term,  while  the  second  term 
corresponds  to  the  collision  term.  The  rate  of  ‘scattering  in'  is  and  the 

rates  of  ‘scattering  out'  is  The  two  rates  have  the  ratio  of 

Boltzmamn  factor.  This  immediately  gives  rise  to  the  equilibrium  distribution 


tt’o(v)  =  e' 


[6.13) 


in  the  absence  of  the  drifting  term  (u  =  0).  When  u  ^  0,  the  Boltzmann  distribution  is  no 
longer  a  zero  mode  of  (6.11).  As  u  deviates  from  zero,  the  pole  of  the  original  zero  mode 
gets  shifted  to  A(u)  <  0.  In  Appendix  D.  we  show  that  A(u)  depends  quadratically  on  u 
as  in  (6.4).  We  found  that  the  diffusion  constant  is  given  by 


_2 _ _ 


(6.14) 


where  d  is  the  dimension  of  the  tight-binding  lattice.  When  the  temperature  is  much  higher 
than  the  band  width  A,  we  may  replace  the  Boltzmann  factors  by  unity.  Thus 

2  <{dev)^>^ 

ft'rfiEvv'^vv'Oev-aev')^'  ■ 

where  the  angular  brackets  denote  average  over  the  energy  band. 
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The  quantity  Ww'  can  be  written  more  explicitly  as 

Wvv'  =(f  r  (6.16) 

*  *  X  1,1' 

where  Fii<(r. x)  was  given  in  (3.19).  For  temperatures  lower  than  the  Debye  temperature, 
we  may  expand  the  exponential  in  the  expression  (3.19)  for  Fn/(r,  x).  The  first  order  term  is 
zero,  due  to  energy  and  momentum  conservation.  A  somewhat  tedious  but  straightforward 
evaluation  of  the  second  order  term  gives  D  where  s  is  the  power  of  with 

which  the  quantity  J(u;)  in  (4.4)  goes  to  zero.  For  short  ranged  interaction,  s  =  5,  so 
D  ~  1/T^,  in  agreement  with  Gogolin[35].  For  a  piezoelectric  interaction,  s  =  3.  so 
D  ~  1/T\ 


6.3  -High  temperature 


What  will  happen  when  the  temperature  gets  high?  Let  us  look  back  to  equation 
(6.10).  The  temperature  dependence  of  £p(u)  and  !3i(u.  A)  are  quite  different.  The  Frank- 
Condon  factor  decreases  exponentially  with  temperature,  so  are  A  and  £p(u).  However. 

the  Frank-Condon  factor  in  Ei(u,  A)  in  (6.7)  will  be  largely  cancelled  by  Fn»(u;,k),  which  | 

i 

grows  exponentijdly  with  temperature  for  I'  =  —  I.  Therefore,  we  can  no  longer  ignore  the  j 
u  dependence  in  Si(u,  A).  As  the  temperature  is  raised  to  a  certain  point,  the  linear  u 
correction  to  Di(0,0)  will  dominate  the  £p(u)  term,  marking  the  transition  out  of  the 
coherent  regime. 

At  very  high  temperatures  we  can  ignore  all  the  exponentially  small  terms,  then  (6.10) 

becomes 

Ac’(v) 


2A 

A2  +u;2 


[-e 


cos(ul)e^^‘*'^^!/>(v  +  k)  ], 


(6.17) 


where  F  =  Fu-  with  1'  =  -1,  and  is  actually  independent  of  1  due  to  the  cubic  symmetry 
of  the  system.  This  equation  has  the  solution  of  (/-’(v)  =contant,  with 

A  =  ,-(^)2y  F(u;,k)— <  cos(ul)  >  (6.1S) 

n  A'^  + 
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where  z  is  the  coordination  number  of  the  lattice,  and  the  angular  brackets  denotes  an 
average  over  1.  Finally,  we  let  A  — *  0+  on  the  right  hand  side  of  the  above  expression,  and 
substitute  (3.19)  into  it,  we  obtain 


A 


“"s'  i 


dr F{t,  0)[—l+  <  cos(ul)  >] 


i6.19) 


where  c  is  the  coordination  number  of  the  tight  binding  lattice,  and  F(r,  x)  is  just  Fii'(  r.  x) 
in  (3.19)  with  1'  =  —1.  The  diffusion  constant  is  then  given  by 


f  / 

exp 

/o  I 

1  ^ 

J(k)  sin^(kli  /2) 
sinh(3^u;k/2) 


COs(ijJicT)  -  1  >  . 


(6.20) 


The  physics  behind  these  arguments  can  be  more  clearly  seen  by  examining  the  self¬ 
energy  operator  (3.20)  in  the  position  representation.  At  high  temperatures,  the  terms 
with  1'  =  —1  dominate,  so  that 


I  u;.k 

g- Jftw/2gik(x-y) 


i6v  v' - ; - ^ — 6 

A  -f-  lijJ 


X  -j"  ILJ 


x.x'  ,y' 


(6.21 


+ 


A  -f  iw 


^x,x'-l-l<5y.y'-H  + 


ik(y-x) 


X  tij 


■*^x.x'  — l^y  .y'  —1 }  • 


It  is  noted  that  only  the  diagonal  transitions  take  place,  i.e.  the  changes  in  x  and  y  are 
the  same.  This  is  in  fact  a  general  feature  of  high  temperature  transport  theory  known  as 
the  qua.si-classicaJ  limit. [36] 
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7.  Conclusion  and  discussion 


In  summary,  we  have  studied  the  question  of  quantum  coherence  for  a  narrow  band 
particle  interacting  with  phonons  and  static  disorders.  We  have  focussed  our  attention 
on  the  transition  probability  P(x/.Xo,f)  for  the  particle  to  move  from  position  Xo  to 
X/  in  time  t.  The  calculation  of  P(x/,Xo,0  and  its  properties  hzis  been  reduced  to  an 
eigenvalue  problem  in  the  phase  space  of  the  particle,  i.e.,  the  space  coordinated  by  the  two 
position  var-ables  appearing  in  the  particle  density  matrix.  The  reduction  is  essentially 
a  renormalization  of  the  tunneling  amplitude  by  the  Frank-Condon  factor,  followed  by 
an  expansion  in  the  small  parameter  A/{h.ujd),  i-e.,  the  ratio  of  the  renormalized  band 
width  to  the  Debye  energy  of  the  substrate  lattice.  The  eigenvalue  problem  is  simple  yet 
general  enough  to  serve  as  a  convenient  basis  for  further  explorations  of  the  various  regimes 
and  their  inter-connections:  low  temperature  versus  high  temperature,  weak  versus  strong 
static  disorder,  two  wells  versus  an  array  of  wells. 

We  have  studied  the  two-well  problem  in  great  detail.  We  have  obtained  complete 
solutions  for  P(x/,Xo,^)  in  the  case  of  low  temperature  but  arbitrary  bias,  and  the  case 
of  zero  bias  but  arbitrary  temperatures.  We  have  demonstrated  how  the  particle  motion 
becomes  incoherent  as  the  temperature  or  bieis  are  increased.  In  the  extreme  cases,  our 
results  agree  with  earlier  ones  in  the  literature. [1], [19]  These  results  provide  useful  guidance 
in  our  study  of  the  many- well  cases  discussed  below. 

For  the  case  where  the  site  energy  disorder  becomes  much  stronger  than  the  tunneling 
amplitude,  we  have  derived  a  meister  equation  for  the  particle  motion.  We  have  found  a 
rather  surprising  result  in  the  regime  defined  by  A  <<  6e  <<  ksT  <<  hujd-  The  diffusion 
constant  is  independent  of  the  disorder  if  the  particle  couples  with  the  phonons  through 
the  piezo-electric  effect,  and  it  increaises  with  disorder  if  the  coupling  is  mediated  by  a 
short  ranged  interaction.  We  also  find  that  the  diffusion  constant  depends  linearly  on 
temperature,  and  depends  quadratically  on  the  renormalized  tunneling  amplitude. 
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In  the  translationally  inv’ariant  case,  we  have  rederived  the  Boltzmann  equation  at 
low  temperatures,  yielding  D  ~  where  s  =  3  for  piezoelectric  coupling  emd  s  =  5 

for  short  ranged  interaction.  At  higher  temperatures,  the  Boltzmann  description  breaks 
down,  zmd  a  transition  to  incoherent  hopping  occurs.  We  have  confirmed  the  Holstein 
formula  for  the  hopping  rate  on  more  rigorous  grounds. 

In  this  paper,  the  term  'incoherence’  is  used  to  indicate  the  Markovian  nature  of  the 
particle  motion  in  positional  basis.  There  is,  however,  a  difference  between  the  incoherence 
occurring  at  high  temperatures  and  that  induced  by  disorder  at  low  temperatures.  The 
physical  principle  operating  in  the  former  case  is  the  quasi-classical  tendency,  the  tendency 
to  move  in  peu'allel  with  the  diagonal  (x  =  y)  of  the  phase  space.  In  the  strongly  disordered 
case,  the  phase  space  diagonal  itself  acquires  a  special  meaning.  It  decouples  from  the 
offdiagonal  sites,  and  becomes  the  main  channel  for  the  particle  motion. 

It  should  be  noticed  that  even  at  zero  temperature,  the  relaxation  rate  T  is  non-zero 
for  the  two-well  case,  bieised  or  unbiased.  One  might  think  that  this  is  the  reeison  for  the 
temperature  independence  of  the  diffusion  constant  reported  in  Ref.26  &27.  In  view  of 
our  results  for  the  case  of  am  infinite  lattice,  one  has  to  conclude,  however,  that  this  is 
not  the  case.  We  expect  that  the  zero-bias  and  zero-temperature  relaxation  rate  should 
go  as  r  ~  J(26^/k)  [see  Equ.(4.29)]  for  a  multi-well  case,  where  6A  is  the  level  spacing  of 
the  renormalized  particle  Hamiltonian.  This  rate  decreases  to  zero  as  the  number  of  wells 
increases  to  infinity.  The  coherence  of  the  particle  motion  makes  it  impossible  to  obtain  a 
finite  diffusion  constant  at  zero  temperature. 

It  is  also  noticed  that  the  coherence  (in  zero  bias)  persists  up  to  a  temperature  of  the 
order  of  the  Debye  temperature  rather  them  the  tunneling  amplitude.  This  is  in  line  with 
the  results  of  Leggett  et  al[19],  who  have  studied  the  general  case  of  a  two- well  particle 
interacting  with  a  heat  bath  through  the  coupling  spectral  function  ~  [see  (4.4)]. 
They  concluded  that  the  transition  to  incoherence  occurs  at  fcgT  =  0  for  0  <  s  <  1.  and  at 
tsT-WsS-j'-forl  <  s  <  2,  where  Wc  is  a  typical  frequency  of  the  heat  bath  modes. 
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It  is  easy  to  see  from  (4.40)  that  the  transition  takes  place  at  ksT  ~  fujjc  for  s  >  2.  For  the 
piezoelectric  coupling  we  have  s  =  3,  and  for  the  short  ranged  interaction  we  have  s  =  b. 
Thus,  a  narrow  band  does  not  necesseirily  imply  an  incoherent  motion.  In  fact,  the  low 
temperature  relaxation  rate  goes  like  T  ~  A*,  whose  ratio  with  the  tunneling  amplitude 
decreases  with  A  for  s  >  1.  indicating  that  the  motion  is  more  coherent  when  the  band 
width  gets  narrower. 

Finally,  it  should  be  pointed  out  that  the  incoherent  hopping  rate  at  high  temperature 
is  proportional  to  the  square  of  the  bare  tunneling  amplitude  A  rather  than  the  renormal¬ 
ized  one,  A.  The  Frank-Condon  factor  cancells  out.  In  Ref.  13,  the  high  temperature  rate 
of  Flynn  and  Stoneham[37]  was  used,  but  the  tunneling  amplitude  was  still  interpreted  as 
the  renormalized  one.  We  believe  this  is  inappropriate. 
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Appendix  A.  Site  energy  change  due  to  displacements  of  lattice  atoms 


Suppose  the  particle  interacts  with  the  lattice  atoms  through  short  ranged  repulsive 
forces.  The  potential  wells  seen  by  the  particle  are  centered  at  the  interstitials  of  the 
atomic  lattice.  Suppose  the  lowest  level  energy  e{L)  of  a  potential  well  depends  on  the 
linear  dimension  L  of  the  well.  A  measure  of  L  mav  be  taicen  ais 


(.-1.1) 


where  Xj  is  the  position  of  a  lattice  atom  bordering  the  well  centered  at  x.  Expanding  to 
first  order  in  the  displacements  X_,  —  X®,  and  reexpressing  them  in  terms  of  the  phonon 
field  operators,  we  then  obtain  the  correction  to  the  site  energy  as  given  in  (2.3).  The 
coupling  function  gilk)  has  the  expression 


'  V 2MuvLl  ‘ 


-X) 


where  is  a  polarization  vector  of  the  phonon  mode.  The  estimate  (2.4)  follows  from  a 
dimensional  analysis  of  the  above  expression. 

Appendix  B.  Reduced  density  matrix  and  influence  functional 

With  a  generalization  of  the  ‘bra-ket’  notations  in  (2.11)  to  the  system  of  our  particle 
plus  phonons,  we  have  in  place  of  (2.12) 


5.1 


where  Ctoi(t)  is  the  Liouvillian  operator  corresponding  to  Htotii)  in  (2.8).  This  is  integrated 
with  the  initial  condition  \ptot{to)  >=  l-^o  >  S|Pp/i(fo)  >.  where  \pph{to)  >  corresponds  to 
the  canonical  distribution  of  phonons  given  in  (2.7),  yielding 


\ptot{i) 


^  ft  fin  ft-. 

>=Y1  / 

n  =  0'tto  Jto  J  to 


dt\ 


n=0“'*o  ''*0 


15.21 


Ch{tn)Ch{tn-x)...Ch[h]\z^>  ^\pph{h)>. 
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where  Cil  corresponds  to  Hr  +  Hl,  and 


iB.3) 


with  everything  else  defined  as  in  the  text.  We  can  then  insert  in  (B.2)  the  complete  sets 
{jcj  ><  Cjl}  for  the  particle  sector,  and  trace  out  the  phonon  degrees  of  freedom,  yielding 

l-t  rtn  f-ti  _ ^ 

<Zf\tT\ptot{t)>=Yl  /  dtn-\...  I  difi  ^ 

„-oJto  Jto  Jto 

X  <*/|£A(iri)|-fn  —  — l|£A(^fi  —  1  )|-^n  — 2 

(BA  I 

n 

X  ex 

j=i 

where  Cii(z)  —<  c|£//,|r  >  now  acts  only  on  the  phonon  degrees  of  freedom.  In  the  above 
and  hereafter,  'tr'  means  tracing  out  of  the  phonon  degrees  of  freedom.  If  the  phonon 
coordinates  are  collectively  denoted  as  R.  then  the  phonon  phase  space  is  coordinated 
as  (R.R').  The  trace  operation  in  the  ‘bra-ket’  notation  is  simply,  trip  >=  'Hr-r'  < 
R.  R'\p>.  Recall  that  the  reduced  density  matrix  is  defined  as 


p{=f,t)  =<:f\tT\ptot{t)>, 


the  expression  (B.4)  is  just  (2.17)  with 

exp(^>[r])  =  tr  >}•  (-5.61 

The  usual  definition  for  the  influence  functional  is,  in  our  language. 


exp('J[r])  =  tr  !pp/,([--].t)>. 


where 


^IPp/.([=Kt)>=  ^/L(-r)|pp/,([c].r)>  . 


The  two  definitions  are  the  same,  if  we  make  the  recognition  of  as  defined  in  (2.16). 
Finally,  the  expression  (2.18)  for  the  influence  functional  has  been  calculated  in  Ref.  15. 
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Appendix  C.  Evaluating  the  influence  functional 


For  a  particular  path  of  n  jumps,  we  can  use  (2.16)  to  carry  out  the  double  time 
integral  in  (2.18).  The  result  is  a  sum  over  contributions  from  the  pieces  of  the  r  -  plane 
partitioned  by  the  times  {tj}.(19]  After  some  rearrangements,  we  have 


$[r]  -  -  J(k)[coth(3ftc.-w/2)Fi  -  iF^] 

k 

n 

j  =  0 

n  ; 

+  cosij,’i,(fj+i  —  fj_r)(F;+i  —  L j)[L' —  L (CM 

j  =  0  r=0 
n 

F2  —  ^  ^  ^  j  (  ^}  + 1  —  ) 

;=:0 

+  ^^Sinu--u(t,  +  l  -  f;-r)(C';+l  -  F,)(V7_,  - 

;=0  r=0 

where  ./(k).  the  U's  and  Vs  are  defined  in  (2.21). 

A  number  of  simplifications  are  in  order.  First,  since  |g(k)|  and  a.'k  are  invariant  under 
k  — +  —  k,  the  first  sum  in  the  expression  for  F\  can  be  replaced  by 

5  E  le,  -  +  \\Uof  +  5IC.V1'.  ,C2i 

“  ;=l 

and  the  first  sum  in  F2  can  be  replaced  by  zero.  These  replacements  do  not  change  the 
value  of  $[r].  Secondly,  since  Ck(tj)  =  0  for  <  0,  we  can  replace  the  lower  limit  to  of  the 
multiple  integrals  over  {t^}  in  (2.17)  by  zero.  Thirdly,  since  co  and  =  -/  are  diagonal. 
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i.e.,  yo  =  xo  and  yn  =  Xn,  we  have  Uq  =  Un  =  0.  We  can  therefore  write 


n— 1 J— 1 


+  +  i  -  tj-r){Uj  +  i  -  Uj)iU"_^  -  L7-r-l  )• 


j  =  l  r=0 


n  —  1  j  —  1 


!C.3'i 


j=l  r=0 


+  sinu-'kft^+i  —  to)(lj j+i  —  L'^;)Vo’' . 


Finally,  we  let  to  — ►  — oo.  This  now  only  affects  the  last  sum  in  F2  in  the  above  expressions, 
whose  contribution  to  <!•[;]  approaches  zero,  because 


does  not  have  sharp  peaics  as  a  function  of  frequency  In  fact,  for  any  finite  time  t. 

are  at  finite  distances  from  Xq.  Therefore,  [Uj+i  —  Uj)Vo  is  a  smooth  function 
of  k  and  vanishes  quadratically  with  k  eis  k  — ►  0.  According  to  discussions  in  Section  1. 
the  quantity  (C.4)  should  behave  smoothly  neeir  a;  =  0.  Finally,  since  J(k)  does  not  have 
sharp  peaks  away  from  k  =  0,  we  expect  the  quatity  (C.4)  to  have  no  sharp  peaks  in 

With  all  these  simplifications,  we  can  now  rewrite  (2.17)  as 


»=  ft  rt„  rt2 

<Zf\p{t)>=y'  /  dtn  /  dtn-\...  /  V 

n=0-^0  '^0  -^0  {z,} 

n 

xexp{A(ro)<i  +  ^  A(r;)(t;+i  -  t^)}  exp($o[-])exp($i[c]), 


where 


1  ” 

=  -  0  E  E  •/(k)coth(.iftu;w/2)i:/;  -  LV,|2. 


j=i  k 


n  — 1 J  —  1 


^  ^  ^  f{^]+  1  ^]-r^  ~j+li  Z}  —  r~  ~)  —  r—l)^ 

j  =  l  r=0 

/(...)  =  -E-^(*‘){coth(JW2)cosu;w{Vi  -tj-r)(Uj+i  -Uj){UJ_,-U;_,_,) 


-  isin^kitj^i  -Uj)[v;_,  -v;_,_,)}. 


(C.6) 


The  factor  exp(<J>[  [-:])  is  just  the  product  of  Fjj>  in  (2.19). 

Finally,  the  time  independent  factor  exp(‘Jo[‘^])  can  be  absorbed  into  the  hopping 
matrix  elements  as 

<:\CkW>  --<z\Ch\z'> 

=  <z\Ch\z'>  exp(-^  Y,  J{^)coth(^hu;i,/2)\U 


Written  explicitly,  we  have 


<x,y|£/,|x',y'  > 

=  i  +  *><y'ly  >  exp(-^  ^  7(k)coth(d;z;.-k/2)|e'‘^"  - 

-  <x|x'><y'|y +  I>exp(-^y'  J(k)coth(i^^’k/2)|e‘‘‘y  - 

2  “  (C.S) 

=  ^  +  l><y'|y>  -  <x|x'><y'ly  +  !>} 

1 

X  exp(-i  V  7(k)coth(dft;*.'k/2)|l  - 

“  k 

where  we  have  used  the  constraints  implied  by  <  x|x'  +  1  >  and  so  forth.  For  an  atomic 
lattice  of  cubic  symmetry,  the  exponential  in  the  last  line  of  the  above  equations  is  inde¬ 
pendent  of  the  direction  of  the  nearest  neighbor  lattice  vector  1,  and  can  be  factored  out 
2is  in  (2.20).  These  arguments  are  still  valid  even  if  the  particle  is  restricted  to  a  subset  of 
the  sites  (potential  wells). 
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Appendix  D.  Evaluating  the  diffusion  pole 

As  u  — +  0,  the  drift  term  can  be  treated  as  a  small  perturbation.  We  first  symmetrize 
the  collision  kernal  by  writing 

0(v)  =  .D.l) 

Then  0(v)  satisfies  the  equation 

A  0(v)  =  (5o0(v)  +  Qi<Z)(v).  iD.2) 


where 

Qo<?(v)  =  ^  VFvv'l^lv')  - 

(£>.3) 

5ev 

Qi<d(v)  =  u— 0(V). 
tn 

We  now  apply  the  Lanczos  iterative  method  to  evaluate  A(u).[38]  We  take  Oolv)  = 
g-J<v/2  starting  state.  We  first  calculate  Qc>o,  project  it  onto  0n,  and  call  the 

remainder  Oi.  We  thus  have 

\  <i>o  =  a  (t>o  0\ ,  ( D A) 

where  by  definition 

a  =  ^  ‘^o(v)0o(v).  (£>.5) 

V  V 

The  coefficient  a  is  actually  zero,  because  (po  is  a  zero  mode  of  Qo  and  the  thermal  average 
of  Qi  =  is  zero.  Next,  we  calculate  (5c>i(v),  project  it  onto  Oq  and  0\,  and  call  the 

remainder  02-  Therefore 

\  <p\  —  b  0Q  A  c  01  A  02>  (Z?-6) 


where 

6  =  <t>o{'v)Q]<Po{'y)/  </>o(v)0o(v), 

V  V 

C  =  ^0o(v)<?iQoQi<^^o{v)/^  0o(v)(5^<^o(v). 

V  V 


{D.7] 


Note  that  while  b  goes  to  zero  quadratically  with  u,  c  approaches  a  finite  constant. 
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If  we  ignore  4>2  in  (D.6),  we  then  have  a  two  dimensional  eigenvalue  problem  (D.4) 


and  (D.6)  to  solve.  The  eigenvalues  are  easily  found  to  be 

Si  ^[c±c(l  +  26/c^)], 


(D.8) 


where  in  the  second  step  we  have  used  the  fact  that  6  <<  as  u  goes  to  zero.  The  solution 
with  the  negative  sign  in  (D.8)  goes  continuously  to  zero  as  u  — f  0,  and  corresponds  to  the 
diffusion  mode.  We  thus  take  A(u)  =  —b/c.  The  expression  (6.4)  and  (6.14)  in  the  text 
then  follows  from  (D.7).  D  is  isotropic  because  of  the  cubic  symmetry  of  the  lattice. 
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FIGURE  CAPTIONS 


Figure  1.  Diagrams  for  the  propagator  in  the  time  domain. 
Figure  2.  The  first  order  diagram  for  the  propagator. 
Figure  3.  The  second  order  diagrams  for  the  propagator. 
Figure  4.  The  self-energy  diagrams. 
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